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Abstract 



This study compares the structure and characteristics of 
tropical synoptic timescale wave disturbances over the western 
Pacific, particularly around the Caroline Islands where six 
upper air stations are present, as contained in radiosonde 
observations and the initial analysis (rOO) and 48-h forecast 
(t 43) of the Navy Operational Global Atmospheric Prediction 
System (NOGAPS) during northern winter of 1991-1992. The 
observed, analyzed and forecast disturbance climatologies are 
all determined from power spectral and cross-spectral analyses 
focusing on the synoptic frequency band encompassing three 
adjacent period points: 4.6-, 5.3- and 6.4-day. The meridional 
wind at 850 hPa is used as the principal parameter for defining 
the disturbances. 

The power spectra of the NOGAPS t 00 analysis agree reasonably 
closely with the radiosonde observations, but the r48 forecast 
shows larger differences. Nevertheless, the observed horizontal 
structure of the disturbances are all well represented in both 
the analysis and the forecast at all levels. These include the 
westward propagation, a 850 hPa wavelength near 33° along 10°N 
and near 38° along the equator, a 200 hPa wavelength near 
50-60°, and a northeast-southwest meridional tilt. In general, 
the analysis and the forecast produce a higher spatial coherency 
than that shown from the radiosonde observations. 

Over the eastern part of the station network the disturbances 
show a strong westward tilt with height similar to that observed 
previously, but at lower levels there is little vertical tilt as 
compared to a significant eastward tilt with height observed 
during some previous summers. Over the western part the 
vertical tilt is eastward throughout the troposphere, which is 
opposite to that of the previously observed summer structure. 
These differences are consistent with the seasonal difference in 
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vertical mean-wind shear due to pronounced monsoon changes. The 
observed zonal variation of the disturbance vertical structure 
is reproduced by the rOO analysis, but the forecast structure 
has a weaker vertical organization over areas where a 
well-organized structure is observed, and a strong organization 
to the west where it is not observed. It appears that the 
model's dynamics produces a downstream propagation that results 
in an artificial wave structure. 

The observed warm core in the upper troposphere is similar to 
previously observed structure during summers and suggests that 
the disturbances are self-supporting energetically. This 
structure is reasonably reproduced by the analysis and forecast. 
On the other hand, the analyzed low-level thermal structure is 
warm core near the surface, opposite to the previously observed 
summer structure. The forecast shows a reverse of this 
structure, from a warm core to a cold core. It is surmised that 
the model's cumulus and planetary boundary layer 
parameter izations may be too strong, causing this reverse. 

Over the northern tropics, the broad scale velocity potential 
shows a southward propagation at 850 hPa and a northward 
propagation at 200 hPa, a pattern consistent with the 
acceleration of the East Asian local Hadley cell observed 
previously during monsoon surges. However, the 48-h forecast 
climatology contains only the 850 hPa propagation. 

Based on the study, a scheme to evaluate the model climatology 
of the disturbances using spectral and cross spectral quantities 
is proposed. 
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1. Introduction 



The analysis and prediction of tropical weather are major 
problems faced by meteorologists due to several reasons, one of 
which is the sparse coverage of radiosonde data, especially over 
tropical oceans. The advances of global NWP models have been 
very useful in alleviating these problems. To make up for the 
lack of observations, the f our-dimensionally assimilated data 
from global models have been used extensively for tropical 
diagnostic studies (e.g., Liebmann and Hendon, 1990; Lau and 
Lau, 1990; Swett, 1993) . The tropical forecasts from the global 
models have also shown useful skill in predicting day-to-day 
weather systems (e.g., Reed et al., 1988; Hogan and Rosmond, 
1992) . Both of these applications of the NWP products are 
extremely important for Navy operations. In the mean time, it 
has been recognized by many investigators that tropical analyses 
and forecasts remain amongst the toughest challenges faced by 
NWP modelers. It is well known that the difficulties in the 
tropical data analysis and initialization, and treatment of 
physical processes such as cumulus convection in the global 
models, continue to cause large errors in the tropics. For 
example, in the early 1980's the tropical divergent circulations 
in all major global operational models were known to be too 
weak. Recent revisions, while resulting in significant 
improvements, still contain large errors, with the ECMWF model 
often having too strong a divergent circulation (Trenberth, 

1992) . 

Many of the error statistics for evaluating tropical analyses 
and forecasts are done in terms of monthly or seasonal mean 
motions where the results describe the distribution of errors of 
mean states. Evaluations of the global model's ability to 
represent tropical synoptic weather systems have been few and 
mostly through case studies by operational forecasters. This 
study is an attempt to evaluate statistically the tropical 
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synoptic disturbances in the NOGAPS analyses and 48-hour 
forecasts over the northwestern Pacific. The region of the 
Caroline Islands was chosen, where several routinely operating 
radiosonde stations allow an evaluation of the analyses and 
forecasts against observations. Spectral analysis was used to 
compare the structures of synoptic disturbances, with 
periodicity around 5 days, as they are deduced from the 
observed, analyzed and forecasted data. The period of study is 
the winter season of 1991-1992. 

2. Data and Methodology 

The data of this study are all retrieved from the FNOC archive 
and may be classified into three categories: 

1) Radiosonde data ( "RAOB" ) : These are reports received 

real-time through the Global Telecommunication System. They are 
the same observations that were used operationally in the 
objective analysis of NOGAPS. Six stations in the Caroline 



Islands and 


vicinity 


are used 


Majuro 


7 . 1 °N 


171. 4°E 


Ponape 


7 . 0 °N 


158. 2°E 


Truk 


7 . 5 °N 


151. 9°E 


Guam 


13 . 6°N 


144. 8°E 


Yap 


9 . 5 °N 


138. 1°E 


Koror 


7 . 3 °N 


134. 5°E 



(These stations are identified in the first figure on a 
geographic map - Fig. 27.) 

2) Analyzed data ("rOO"): The objectively analyzed fields 

in NOGAPS. The data are on the 2.5° x 2.5° latitude-longitude 
grid. 
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3) Forecast data ("rhh"): The NOGAPS forecasted fields at 

hh hours. For most of the study the 48-h forecast (t 48) is 
used. The data are on the 2.5° x 2.5° latitude-longitude grid. 

All data are twice daily, from 12Z 11 November 1991 to 00Z 17 
March 1992, a total of 128 days or 256 time points. Missing 
data are filled by linear interpolation of the nearest values. 
Time series of radiosonde with coverage less than 80% are not 
used . 

Spectral and cross-spectral analysis are used to determine the 
characteristics and structure of synoptic timescale fluctuations 
contained in the data, as well as to compare the NOGAPS analyses 
and forecasts against radiosonde observations. The spectral 
analysis is done by partitioning each demeaned and detrended 
time series into seven ensembles with a cosine taper of 10% 
applied to each end before fast Fourier transformation. The 
resultant frequency resolution is 0.03125 cycles/day (32-day 
period) , and the Nyquist frequency is 1 cycle/day. The (32 
day) -1 point is the lowest frequency resolved, and because of 
the detrending, a spectral peak at this frequency is considered 
to represent all variances beyond the monthly timescale, to 
include those of intraseasonal and possibly seasonal 
oscillations . 

The determination of the disturbance structure is similar to the 
approach used in previous spectral studies on the easterly waves 
of the tropical western Pacific using radiosonde data (e.g., 
Wallace and Chang, 1969; Chang et al., 1970a; 1977b). 

Inspection of all spectral distributions revealed that 
significant variances often occur at seasonal, intraseasonal, 
quasi-biweekly and synoptic timescales. Cross-spectral analysis 
is then used to focus on the synoptic scale oscillations, whose 
variances are mostly concentrated in a frequency window centered 
at period 5.3 days. The cross spectra, in terms of coherence 
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squares and phase differences, of the three adjacent frequencies 
corresponding to 4.6, 5.3 and 6.4 day periods, are then averaged 
to represent the cross relationship between various pairs of 
variables of the synoptic scale. This frequency band will be 
called the synoptic band or the "5-day" band. The cross spectra 
are used to infer the horizontal, vertical and the dynamic and 
thermodynamic structures of the 5 -day disturbances. 

The degrees of freedom for the 5-day band is estimated to be 42, 
and the corresponding coherence squared threshold for the 95% 
significance level is estimated to be 0.14. Phase difference 
values with coherence squared above 0.14 will be considered 
significant and are used for inferring the disturbance 
structure. (When an individual frequency is considered, the 
degrees of freedom is about 14 and the 95% significance level is 
at coherence square > 0.43.) 

The meridional wind oscillations are generally a better 
indicator of tropical synoptic-scale disturbance activities 
except at the equator. This is because, for the duration of a 
synoptic time period, the tropical winds are normally dominated 
by the zonal component with a high constancy, and the meridional 
wind tends to oscillate around a very small time mean. Based on 
previous spectral studies of westward propagating synoptic 
disturbances in the tropical western Pacific using radiosonde 
data, the meridional velocity at 850 hPa (V 35 Q) is used as the 
most important variable in detecting the synoptic disturbances. 
850 hPa is the level where the structure of these disturbances 
has been found to be most clearly defined, and it is also the 
level of maximum data coverage in the lower troposphere because 
of the existence of many pilot balloon observations and cloud 
base satellite winds. In many cases vg 5 Q is used as the base 
time series to cross with other series in order to determine the 
structure of disturbances. In the upper troposphere, the 
meridional velocity at 200 hPa, where maximum amount of 
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upper-level satellite cloud winds and commercial aircraft wind 
data are located, is considered the base series. 

The area of study covers the entire tropical western Pacific. 

All horizontal maps cover the region from 100°E-180°E and 
10°S-30°N. Because of the wide longitudinal span of the region, 
the term "central-western Pacific" will occasionally be used to 
refer to the eastern half of the diagram (140°E-180°E) . While 
the main focus of the study is the synoptic westward propagating 
disturbances in this region, the South China Sea is included in 
the diagrams because of its importance in the winter monsoon 
cold surge events that are known to dominate the synoptic scale 
fluctuations in the lower troposphere over the western Pacific. 
For brevity, the term "easterly waves," where "easterly" refers 
to propagation from the east, may be used occasionally to refer 
to the tropical westward propagating synoptic disturbances. 

3. Comparison of power spectral characteristics 

The power spectra of the vg 5 Q time series for all three 
categories of data are calculated first to compare the spectral 
characteristics of the NOGAPS analysis and 48-hour forecasts 
with those of the radaosonde observations. Overall, prominent 
oscillations (with large spectral power or power peak) usually 
belong to one of the following three timescales: the 
seasonal/intraseasonal, the 10-15 day period, and the synoptic 
(centered at 5-day) period. We will focus on the synoptic 
period in this study. 

a. Meridional winds at 850 hPa 



Figs. 1-6 show the comparisons of the V 350 power spectra between 
RAOB, TOO analysis and r48 forecast at each of the six stations, 
arranged from the easternmost station (Majuro) to the 
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westernmost station (Koror) . The TOO and t48 series are based 
on the grid point closest to the station. In each figure, the 
top panel (a) displays the power distributions of the RAOB 
(heavy solid), TOO (thin solid) and r48 (dashed) series, with 
all curves normalized with respect to the maximum value within 
the diagram. The middle panel (b) plots the coherence squares 
and phase difference between the RAOB and rOO series, while the 
lower panel (c) displays the coherence squares and phase 
difference between the RAOB and r48 series. We will discuss the 
tOO-RAOB comparisons in all six figures first, followed by the 
T48-RAOB comparisons. 

The spectral distribution in the RAOB ' s varies between stations, 
but in general two local spectra peaks may be identified. At 
the low frequency end, there is a peak near 32 days. Since the 
linear trend has been removed from the time series, this peak 
may represent both the seasonal trend and the intraseasonal 
oscillations, including possibly the Madden-Julian Oscillation. 
This low-frequency peak appears in all six stations except 
Majuro, where significant power is nevertheless indicated. At 
the synoptic frequency range, a local peak near five days is 
indicated in the RAOB data of all six stations, with Ponape, 

Guam and Koror showing substantial variance in the distribution. 
Between these two frequencies another power peak is often 
indicated, but the frequency location and the relative variances 
vary considerably from one station to another. At Majuro the 
RAOB's have a dominant peak near the 10-day period, 
significantly larger than the power near 5-day and overshadowing 
that of the the 32-day. At Truk, Guam and Yap, an indication of 
a peak is found near 8-day. At Ponape and Koror, the 
lower-frequency peak covers a broad band from 10- to 32-day. 

i) Vg 50 / TOO vs RAOB 



9 



The shape of the TOO power spectra is in reasonably good 
agreement with the RAOB spectra at Ponape, Truk, Guam and Koror. 
At Ponape and Truk the NOGAPS analysis overestimates slightly 
the 5-day variance and underestimates slightly the lower 
frequency variances. At Guam the 5-day TOO variance is close to 
that of the RAOB's, and at Koror the analysis underestimates 
slightly the observed 5-day variance. The coherence squares at 
all four stations are mostly very high and the phase differences 
almost uniformly small throughout the synoptic and lower 
frequency range, indicating a good objective analysis 
climatology of the synoptic scale variabilities. 

Large differences occur at the two other stations. At Majuro 
the TOO spectral distribution is quite different from the 
RAOB's, especially at the lower frequencies. Near the 5-day 
band the variance magnitudes are comparable but the shape is 
different, with the RAOB peak at 4.6-day but the analysis peak 
at 6.4-day. The coherence values are highest in the 4- to 5-day 
band, indicating that the TOO analysis is still reliable for 
representing the synoptic scale fluctuations. At other 
frequencies the coherences are lower than the four 
aforementioned stations, but they are all substantially above 
the significant threshold. The phase difference between the 
analysis and RAOB data is very small near 5-day, but 
non-negligible differences can be seen in adjacent frequencies. 
Near the 8-day point, TOO leads the RAOB oscillations by about 
45°, indicating the analyzed disturbances move through this 
station at a speed 1/8 times faster than the observed. Near the 
3-day point TOO lags RAOB slightly. These two points are also 
where the coherence squares fall off substantially. 

At Yap the magnitude of the entire TOO spectrum is only a 
fraction of the observed. This is apparently because of the 
large RAOB variances with its maximum value nearly seven times 
(maximum Vg 5 Q fluctuation amplitude nearly 2.5 times) that of 
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Koror, a station about five degrees to its southwest. This 
large fluctuation seemed to be ignored by NOGAPS as local 
anomalies, and the model chose to analyze towards the large 
scale flow, suppressing the RAOB fluctuations at Yap. 
Nevertheless, the two largest rOO peaks, at 32-day and 4.6-day, 
are in agreement with the RAOB. These are also the only 
frequency points where both the phase difference is small and 
the coherence squares is > 0.4. 

ii) V 350 / t48 vs RAOB 

Figs. 1-6 show that in all cases the correspondence between the 
r48 and RAOB V 350 spectra is less than that between rOO and 
RAOB. At Majuro (Fig. la) the t 48 synoptic peak occurs at 
6.4-day rather than the observed 4.6-day. This is the same peak 
frequency as shown in the rOO spectrum, but unlike the rOO data 
where the magnitude of the 4-5 day variance is comparable to 
that of the RAOB, the r48 variance is less. The coherence 
squares (Fig. lc) are mostly lower also, with larger phase 
differences. In fact, in the rOO-RAOB coherence distribution 
the values at the observed 4-5 day peak are close to the maximum 
in the entire frequency domain, but for T48-RAOB the lower 
coherence values at this period band is relatively low within 
the frequency distribution. The maximum coherence between t 48 
and RAOB occurs at 6.4-day, which is the fictitious peak present 
in the NOGAPS analysis and forecast, but not in the radiosonde 
data . 

At Ponape (Fig. 2) the good correspondence between rOO and RAOB 
spectral distributions is not reproduced by T48. An increase in 
the power between 5- and 10-day occurs in the r48 spectrum, and 
the observed 4.6-day peak is shifted both directions to 4-day 
and 6 . 4 -day. At Truk (Fig. 3) the observed 5. 3 -day peak is 
shifted to 6.4-day. At both Ponape and Truk the observed 
low-f requency power is substantially reduced. Thus a 
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substantial part of the observed seasonal trend is lost in the 
t 48 forecast. The T48-RAOB cross spectra at both stations also 
indicate slightly lower coherence squares and slightly larger 
phase differences compared to the rOO-RAOB cross spectra, 
although their values are still reasonable enough to indicate a 
fairly good agreement in the lower frequency and synoptic 
oscillations. 

At Guam (Fig. 4) the power of the r48 spectrum around the 10-day 
frequency band is significantly larger than observed. Otherwise 
the spectral distribution is similar to the rOO spectrum. The 
good agreement of the synoptic peak around 5-day between rOO and 
RAOB is changed slightly, with a 20% decrease of power and a 
shift of the peak from the observed 5.2-day to 4.6-day. There 
is a reduction in coherence squares but the values remain 
significant throughout the lower- to synoptic frequency domain. 
Particularly, the coherence square around 5-day is very high, 
indicating very good agreement for the 5-day oscillations. The 
phase differences are again only slightly larger than the near 
zero rOO-RAOB values. 

As is the case with the TOO data at Yap (Fig. 5) , the t 48 
oscillations have much smaller power compared to the very large 
RAOB fluctuations, and the coherence and phase data indicate 
almost no agreement. On the other hand, the t 48 data show a 
predominant spectral peak at 10.7-day that is not present in the 
rOO analysis, where the peak power occurs at the lowest 
frequency (32-day) . 

The deterioration of the t 48 spectrum from rOO at Koror (Fig. 6) 
is similar to that at Ponape, except it is somewhat more severe. 
The good agreement of the 5.2-day peak between TOO and RAOB 
spectra is totally lost, and is replaced by a forecasted 
10.7-day peak. In the RAOB data, 10.7-day is not a peak but is 
on the edge of a broad low-frequency peak, with substantial 
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power. The coherence between T48 and RAOB is generally much 
reduced except around 11-16 days, and the phase difference is 
small only for these and lower frequencies. Since Koror and Yap 
are closely located to each other, it is possible that the 
10.7-day peak present in Yap's r48 forecast is related to the 
one forecasted by Koror. 



b. Zonal winds at 850 hPa 



Comparisons of the RAOB, TOO and t 48 ug 5 Q spectra for the six 
stations are shown in Figs. 7-12. Again, the comparisons 
between RAOB and TOO will be discussed first, followed by the 
comparisons between RAOB and t48. 

The observed (RAOB) Ug 5 0 power distributions are typically 
different from those of the corresponding Vg 5 0 , and the 
magnitude of the Ug 5 Q variances are generally larger than the 
v 8 50 variances except at Majuro and Yap where the two are 
comparable. Table 1 lists the ug 5 Q and vg 5 Q spectral peaks in 
the lower-frequency band and the synoptic band (close to 5 days) 
observed in the RAOB data: 

Table 1 

v 850 u 850 



Station 


low-f rea 


synoptic 


low-f rea 


synoptic 


Majuro 


10.7 


4 . 6 


>32 


weak 


Ponape 


16 


4 . 6 


>32 


6.4 


Truk 


>32 


5.3 


>32 


4.2 


Guam 


>32 


5.3 


>32 


4 . 6 


Yap 


>32 


4 . 6 


16 


weak 


Koror 


16 


5.3 


16 


weak 


Compared to the 


observed v 


spectra. 


the power in the ' 



is more concentrated in the lower frequencies. In the synoptic 
band, the peaks for Ug 5 0 and vg 5 Q at each station are usually 
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different, and a relatively higher fraction of power is found in 
the v spectra than in the u spectra. This is consistent with 
previous findings (e.g., Wallace and Chang, 1969), and is an 
indication that the lower-frequency seasonal or intraseasonal 
oscillations are best represented by the u winds, and the 
synoptic timescale oscillations are best represented by the v 
winds. We will therefore use the v component as our primary 
time series to deduce the structure of synoptic scale 
disturbances in the RAOB , analysis and forecast data. 

i) Us50/ 700 vs RAOB 

In general, at most stations the agreements between the spectra 
of rOO us5o and its corresponding RAOB spectra (Figs. 6-12) are 
better than those of rOO vgso (Figs. 1-6) , particularly the 
western stations of Truk, Guam, Yap and Koror. At these 
stations the spectra of rOO ug 5 Q are almost a reproduction of 
the respective RAOB spectra, which is not the case for the TOO 
Vg5C) spectra. The phase differences between TOO and RAOB at 
these stations are near zero throughout the lower-frequency and 
the synoptic band, and are marked by very high coherence values. 
This is true even at Yap (Fig. 11) , where, in contrast to the 
rOO Vg 5 Q whose variance is only a fraction of the observed large 
fluctuations and the coherence and phase values show poor 
agreement with RAOB (Fig. 5) , the TOO ug 5 Q is in good agreement 
with the observed ug 5 Q in spectral variance and distribution, 
and in coherence and phase values. 

At the eastern end of the station network, Majuro (Fig. 7) shows 
the largest difference between TOO ug 5 Q and RAOB, but the 
general agreement is still somewhat better than the 
corresponding chart for vg 5 Q (Fig.l). This is particularly true 
at the lower frequencies where the coherence square and phase 
difference are better also. The low-frequency spectral peak at 
Majuro occurs at 32-day for both RAOB and the TOO analysis, 
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indicating significant seasonal and/or intraseasonal power. 
However, the rOO peak is rather confined to the lowest frequency 
and does not indicate the significant spectral power at 16-day 
that is clearly seen in the RAOB. The synoptic band at Majuro 
is different: both TOO v and u at 850 hPa show peaks at a 
frequency different from the observed, but the magnitude of rOO 
v 850 power is comparable to that of RAOB while the rOO Ugso 
power is less. The coherence and phase values of the V 350 
t OO-RAOB also show better agreement than the corresponding ugso 
values . 

At Ponape the rOO Ugso spectral shape (Fig. 8 ) resembles RAOB 
better than the rOO vg 5 Q (Fig. 2) , especially at the lower 
frequency where the coherence square with RAOB for ug 5 Q is very 
high. However, similar to Majuro, the Ponape rOO spectrum fails 
to show the observed significant power at 16-day that is part of 
the broad low-frequency peak. At the synoptic band both TOO 
ug 5 Q and vg 5 Q have the correct frequency for a spectral peak 
that is the same as the corresponding RAOB peak, with the former 
underestimating the observed variance, while the latter 
overestimates. The TOO ug 5 Q-RAOB coherence square is also very 
high, showing a value close to 1 at the dominant synoptic peak 
of 6.4-day. 

ii) Ug 5 Q, r48 vs RAOB 

As expected, the comparison of spectra between r48 ug 5 Q and the 
corresponding RAOB data in Fig. 7-12 show that there is some 
degradation from TOO, but the extent is generally small, and is 
less than the degradation of the vg 5 Q t48 forecasts. At Majuro 
(Fig. 7) the t48 Ug 5 Q spectrum shows the same lowest frequency 
power at 32-day as rOO, but the variance is much higher than 
RAOB. The coherence is also much reduced. The entire t48 power 
spectrum appears similar to the RAOB spectrum. It now produces 
a peak with comparable magnitude to the RAOB peak at 8 -day, 
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seemingly an improvement over the r 00 spectrum. However, the 
coherence square with RAOB is very low, indicating that the 
relationship between the t48 Uqso forecasts and RAOB at Majuro 
is rather poor. 

At Ponape (Fig. 8) the low-frequency power distribution is 
similar to rOO with power concentrated at the lowest frequency 
(32-day). Unlike the t 00 data, the observed 6.4-day peak 
appears to be shifted to 8-day in the r48 spectrum. The 
coherence squares with RAOB for the band with frequency lower 
than the 5-day period are around 0.5, about half of the near 
perfect values given by the rOO data but still significant. The 
phase relationship shows small but discernible differences with 
RAOB, which is not the case for the rOO data. 

At the western stations the r48 Uqso forecasts generally 
correspond fairly close to the rOO and the RAOB spectra. The 
exceptions are: at Truk the observed 6.4-day peak is shifted to 
8-day in 748, a feature similar to the neighboring station of 
Ponape; at Guam and Yap the low frequency variance is somewhat 
overforecasted; and at Koror the 16-day spectral peak is shifted 
to the lowest frequency of 32-day. In general the coherence 
squares with RAOB in the low-frequency and synoptic band is high 
or at least significant. The phase differences are similar to 
Ponape, small but more discernible than those shown by the TOO 
analyzed data. 

c. Meridional winds at 200 hPa 

The comparison between the RAOB and the r00 and r48 V 200 spectra 
for the stations are shown in Figs. 13-18. At all stations a 
dominant low-frequency spectral peak occurs at the 32-day 
period, which may or may not coincide with those observed at 850 
hPa . The frequency locations of the synoptic peaks at the two 
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levels are more in agreement, although the magnitude of the 
relative and absolute variances may be different. 

i) V200/ 700 vs RAOB 

At Majuro the observed 200 hPa peaks are at 32-day, 8-day and 
4.6-day (Fig. 13). The 4.6-day peak is the same frequency as 
that found at 850 hPa, but the major 850 hPa peak of 10.7-day is 
missing at 200 hPa and is replaced by a power minimum. The TOO 
power spectrum simulates the 8-day peak very well and is 
supported by high coherence and near perfect in-phase 
relationship with RAOB. At the other two peak frequencies the 
t 00 analysis underestimates the variances, and the lower 
coherence squares (and the small but noticeable phase difference 
at the lower frequency) indicate that the correspondence with 
RAOB are not good. 

At Ponape (Fig. 14) , the dominant peak occurs at the lowest 
frequency (32-day period), and relative to it all other peaks 
appear to be rather minor. The major Vg 5 Q peak at 4.6-day is 
one of the minor ones and is not simulated well by the TOO 
analysis. The coherence and phase values between r00 and RAOB 
can only be considered as indicating a moderate correspondence 
between the two time series at this frequency. At Truk (Fig. 

15) an 8-day peak is observed at 200 hPa and seems to correspond 
with one of the Vg 5 Q spectral peaks, but the other synoptic Vg 5 Q 
peak at 5.2 day is not observed. The TOO spectrum only 
simulates the observed 32-day peak, although the coherence 
square around the 8-day period is quite high. The phase 
difference throughout the synoptic and low-frequency bands is 
reasonably small. 

At Guam, Yap and Koror (Figs. 16-18) two observed peaks occur: 
the commonly-seen 200 hPa peak at 32-day and another at 
10.7-day. Both are detectable in the TOO analysis at all three 
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stations, although the signal for the 10.7-day peak is rather 
weak at Yap and Koror. The RAOB power spectrum over the 
synoptic band is fairly closely simulated by tOO also, although 
the minor peaks around 5-day are not exactly reproduced at Guam 
and Yap. For these two stations the coherence squares are very 
high at the 32-day and 10.7-day peaks, and are quite high 
elsewhere. At Koror the coherence square at 10.7-day, although 
still significant, is a minimum. The phase differences are 
reasonably small throughout the low and synoptic frequencies. 

ii) v 200/ t48 vs RAOB 

Compared to the TOO analysis, the r48 forecasts shown in Figs. 
13-18 generally show a degradation in the loss of power and a 
deterioration of the coherence and phase values. The largest 
loss of power occurs in the eastern stations of Majuro (Fig. 13) 
and Ponape (Fig. 14). The only observed feature that is 
generally reproduced is the 32-day peak. At Majuro the r48 V 200 
peak is too confined to the lowest frequency, and at Ponape it 
is substantially underestimated. At other stations the lowest 
frequency peak is rather well simulated by the r48 spectra. 
However, the good simulation of the spectrum does not guarantee 
a good cross-spectrum result. This is clearly seen at Yap (Fig. 
17) , where the 32-day coherence square between r48 and RAOB V 200 
is rather low even though the two spectral peaks at 32-day are 
almost identical. With the exception of Guam, the phase 
difference between t48 and RAOB time series at all other 
stations are large over the low-frequency and synoptic band. 

d. Zonal winds at 200 hPa 
i) U 20 O/ r0 ° vs RAOB 

The observed spectra U 200 at a H stations (Figs. 19-24) shows 
power concentrating in the lowest frequency, and is well 
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represented by the rOO spectra at all stations with nearly 
identical powers, very high coherence squares and small phase 
differences in the rOO-RAOB cross-spectra. There is practically 
no spectral peak near the 5-day band, except a very slight 
indication at Truk. On the other hand, Guam and Yap (Figs. 22 
and 23, respectively) show a secondary peak at 8-day. Only the 
Guam 8-day peak, which appears more significant than the one at 
Yap, is simulated by the rOO spectrum with an underestimate of 
the observed variance by about 1/3. 

ii) U 20 O/ t48 vs 

Except at Guam, the r 48 power spectral distributions in Figs. 
19-24 generally resemble the observed U 200 spectra with power 
concentrated in the lowest frequency. However, the excellent 
reproduction of the magnitude of the variance by the rOO 
analysis no longer exists. At almost all stations the variance 
represented by r48 U 200 i s substantially reduced, often by 50% 
or more. Only at Koror (Fig. 24) is the reduction less severe, 
which is probably a result of the overestimate of the variance 
in the rOO data. At Guam (Fig. 22) the r48 spectral shape is 
different, with a broad peak occupying the band between 8-32 
day, which is not observed or analyzed. In general the T48-RAOB 
coherence and phase values deteriorate significantly from the 
comparable r 00 values. 



4. Zonal structure of the synoptic wave disturbances 
a. Longitudinal cross spectra of V g 5 Q along 7.5°N latitude 

Eecause the v 8 5 0 power spectra generally show the clearest 
signals of the synoptic (5-day) band disturbances, the zonal 
structure for these disturbances in terms of wavelength and 
phase velocity is determined by the cross spectra of V 350 along 
the same latitude. The Caroline Islands stations are more or 
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less lined up longitudinally along the 7.5°N zone with Truk in 
the center. We therefore use Truk V 350 as the base series to 
perforin cross-spectral analysis with vg 5 Q at other stations to 
determine the observed zonal structure in this latitudinal zone. 
This is done at three frequencies corresponding to the periods 
of 4.6, 5.3 and 6.4 days. The results are shown in Fig. 25. In 
each of the panels the phase difference of V 350 at a station 
with respect to Truk is plotted by the leading letter of the 
station name, and the coherence square is indicated by a star 
sign at the same longitude. In general the coherence falls out 
as the distance of a station from Truk increases, but the values 
all meet the 95% significance level test. The phase difference 
shows consistently a westward (easterly) propagation with a 
wavelength near twice the distance between Majuro and Ponape, or 
about 3300 km. This estimated wavelength increases with the 
period. It is somewhat shorter at 4.6-day and somewhat longer 
at 6.4-day, to give a easterly phase speed that is fairly 
steady. It may be interesting to note that this 
semi-nondispersive characteristic has been noted in 
longitude-time series of satellite cloud data over the tropical 
western Pacific (Chang and Lim, 1982) . 

The signals of the synoptic wave disturbances are most prominent 
in the vg 5 Q time series. Because of the vague signals in other 
parameters, the station RAOB data will not be used as the basis 
to determine the zonal structure. For these parameters the 
comparison will be restricted to the general characteristics 
obtained from the TOO and t 48 data. 

i) Comparison of TOO vs. RAOB 

Also plotted in Fig. 25 are the coherence squares and phase 
differences calculated from TOO V 350 at grid points along 7.5°N, 
with TOO V 350 at the point nearest Truk as the base series. It 
can be seen that the phase line of the TOO cross-spectra agrees 
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quite well with the phase points calculated from the station 
RAOB data at all three frequencies, particularly to the east of 
Truk. To the west, data from the grid point closest to the Guam 
longitude agrees well with RAOB, even though Guam is the 
farthest station from 7.5°N. Apparently the wave disturbance 
has sufficient meridional scale without a steep tilt to allow 
this result. At 4.6-day the grid point data near Yap has a 
smaller phase difference (implying shorter wavelength) than that 
indicated by the station RAOB, and at 6.4-day the grid point 
data near Koror has a larger phase difference (implying longer 
wavelength) than the RAOB data. Thus the r 00 analysis minimizes 
the increase of the estimated wavelength as the wave period 
increases. Namely, r 00 gives a more constant estimate of the 
zonal wavelength within the synoptic band than the station data. 
The TOO coherence squares also show a consistent fall off away 
from the base point near Truk. They are in general higher than 
the RAOB values, which can be expected because of the spatial 
objective analysis. 

ii) Comparison of r48 vs. RAOB 

Fig. 26 is similar to Fig. 25 except here the cross spectra are 
determined from the t48 time series. The agreement between 
these results and the station RAOB or r00 analysis time series 
are fairly good, except for regions west of Yap, where for both 
the 5.3- and 6.4-day periods the consistent westward propagation 
is lost in the t48 fields and the coherence squares are 
meaninglessly low. An examination of Figs. 11-12 indicate that 
the r48 forecast of vg 5 Q is very poor near Yap and Koror. 



b. Longitudinal cross spectra of v qka as a function of latitude 

The good agreement between TOO and RAOB vg5Q longitudinal 
structure suggests that an estimate of the phase propagation of 
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the synoptic wave disturbances based on rOO vg 5 Q is worthwhile, 
and that a comparison between t48 and rOO V350 may provide a 
meaningful evaluation of the models performance of forecasting 
the propagation characteristics of the tropical synoptic waves 
in the vicinity of the Caroline Islands. In this section the 
zonal structure at each latitude will be calculated from the 
phase differences and coherence squares of the grid point data, 
for both rOO and t48, with the 150°E data as the base series. 

The 150°E base longitude is chosen because it is close to the 
longitude of Truk, the base station used in the RAOB 
comparisons. These cross spectra will be referred to as the 
inter-x cross spectra, and are computed for different latitudes 
and presented on a two-dimensional horizontal map. 

i) Inter-x rOO vg 5 Q 

Fig. 27a shows the inter-x cross spectra for rOO. Here at each 
grid point the cross spectra between the vg 5 Q time series and 
that at 150°E of the same latitude is plotted as a vector, whose 
direction indicates the phase difference and whose magnitude is 
proportional to the coherence square. For coherence squares 
less than the 95% significance level, the vector is not plotted. 
Zero phase difference is indicated by a vector pointing 
northward (0°), and as the vector rotates clockwise the positive 
phase difference increases. For vectors pointing to the east, 
the base series (vg 5 Q at 150°E) lags, and for vectors pointing 
to the west the base series leads. Depending on the position of 
a grid point relative to the base longitude, the direction of 
propagation can be determined, and a zonal wavelength can be 
inferred from the phase differences. 

In Fig. 27a the phase propagation characteristics of the 
synoptic band waves as a function of latitude is clearly 
depicted. At and north of 17.5°N eastward propagation can be 
found over a large longitudinal span, with a wavelength of about 
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60° longitude. This is apparently the signals generated by the 
winter monsoon synoptic systems (cyclones and anticyclones) as 
they move out of the Asian continent into the western Pacific 
(eg. Boyle and Chen, 1987) . At and south of 15°N the easterly 
wave disturbances are fairly well defined by the westward phase 
propagation, although significant coherence values are limited 
to a more confined area within the tropical western Pacific east 
of 135°-140°E. The implied zonal wavelenth also varies somewhat 
with latitude. The phase vectors to the east of the base 
longitude indicate a half-wavelength that changes from about 15° 
longitude at 10°N to about 17.5° longitude at the equator. The 
phase vectors to the west of the base longitude suggest a 
slightly longer half-wavelength, but the coherence values fade 
below the 95% significance level. Overall, the phase structure 
suggest a zonal wavelength of about 33° longitude at 10°N, in 
agreement with the results shown in the tOO-RAOB inter-x 
comparison (Fig. 25) . Along the equator the wavelenth is about 
38° longitude. Immediately to the south of the equator the 
wavelength appears to be even longer, especially to the west of 
the 150 °E base longitude, however the very low coherence 
prevents an estimate of the wavelength there. 

ii) Inter-X r48 vg 5 Q 

Fig. 27b shows the inter-x cross spectra for r48 vg 5 Q, again 
with 150 °E as the base longitude. The result is an excellent 
reproduction of the corresponding diagram of TOO (Fig. 27a). 

The entire latitudinal distribution of horizontal wavelength and 
coherences are very well represented by the 48-h forecast data, 
even including some of the details that are somewhat remote from 
the base longitude. For example, the cross spectral vectors to 
the south of Sumatra in the lower-left corner of the diagram are 
nearly identical to the results obtained in the TOO diagram. 
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There are a few above-significance level phase vectors in the 
r48 vg 5 o results that do not show up in the rOO results. These 
are in the 130-140°E, 10°S-equator region southwest of New 
Guinea and indicate a westward propagation of rather long 
wavelength. (A half-wavelength of near 40° longitude can be 
inferred from the phase vectors) . This long wavelength is in 
the range of the mixed Rossby-gravity wave that was observed in 
the equatorial upper-troposphere (Yanai and Maruyama, 196x) . It 
is worth noting that this indication of the planetary-scale 
waves is undetectable in the r 00 analysis (Fig. 27a) . 

c. Longitudinal cross spectra of u g^n 
i) Inter-X TOO ug 5 Q 

Fig. 28a shows the inter-x cross spectra for the zonal wind at 
850 hPa with 150°E as the base longitude. In the northern 
domain the eastward propagation characteristics, including the 
coherence and phase values and their inferred wavelength, are 
very similar to the meridional wind results (Fig. 27a) . 

However, the domain of eastward propagation here covers all 
latitudes north of 12.5 °N, five degrees to the south of the 
lowest latitude of eastward propagation determined from the vggQ 
inter-x cross spectra. Thus it appears that the signals of 
tropical easterly wave disturbances are mainly in the v 
component, and are masked by the subtropical westerly 
propagating systems whose u fluctuations remain strong in the 
12 . 5 ° — 17 . 5 °N zone. 

In the equatorial latitudes the u results are more different 
than the v results found in Fig. 27a. There is an indication of 
westward propagation, but the coherence squares are 
significantly lower between the equator and 10°N, except in the 
immediate vicinity and slightly to the east of the base 
longitude. Based on this partial wave structure, the implied 
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wavelength should be longer than that shown in the meridional 
wind results. At the equator and southward, the results suggest 
an even longer wavelength in the westward propagation over 
Indonesia. The wavelength is about 120° longitude, or 
wavenumber three. This clear planetary-scale wave propagation 
in the u cross-spectra contrasts sharply with the lack of a 
coherent propagation signal in the v cross spectra in the 
equatorial zone. 

ii) Inter-X T48 ug 5 Q 

The r48 Ug 5 Q inter-x cross spectra shown in Fig. 28b resemble 
closely the TOO results (Fig. 28a) , with eastward propagating 
synoptic-scale waves north of 12.5°N, and partial structure of 
westward propagating wave south of 10 °N. The near wavenumber 
three planetary wave is again found in the latitudes south of 
the equator. 



d. Longitudinal cross spectra of v opp 
i) Inter-X TOO V 200 

At 200 hPa the TOO V 200 inter-x cross spectra with 150°E as the 
base longitude again displays the eastward propagation in the 
northern latitudes and westward propagation at and south of 10 °N 
(Fig. 29a) . However, the westward propagating wavelength is 
close to 50° longitude to the east of 150°E, and somewhat longer 
and less well defined to the west. Thus, the upper tropospheric 
wavelength is longer than that at 850 hPa. A longer upper 
tropospheric wavelength has been reported by many previous 
observations based on radiosonde data (e.g., Chang et al, 1970; 
Chang and Miller, 1977), and suggests that the vertical wave 
structure may be complex. The lower and upper tropospheric 
westward propagating disturbances may be far from tightly 
coupled. 
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ii) Inter-X r 48 V 200 



Fig. 29b shows the inter-x cross spectra for t 48 V2oo* In 
general, it is in good agreement with the patterns of rOO (Fig. 
29a) . The difference occurs mainly east of the 150°E base 
longitude, where between 10°N-20°N the eastward propagating 
waves show better organization than the TOO results with higher 
coherence values. South of 10 °N the westward propagating waves 
are less organized than the TOO results, having lower 
coherences . 

e. Longitudinal cross spectra of u onn 

i) Inter-X TOO U 200 

Fig. 30a shows that the rOO analyzed U 200 again indicates 
eastward propagation in the north and westward propagation in 
the south. However, the tropical westward propagation, which is 
evident at 15 °N and southward, has a much longer wavelength than 
that indicated by the TOO V 200 inter-x cross spectra (Fig. 29a) . 
Based on the few significant phase vectors over Borneo, Sumatra 
and Java (in the lower-left corner of the diagram), the 
quarter-cycle phase difference may be placed near 120°E, which 
gives an estimate of a wavelength of 120° longitude. This is 
the same estimate as the 850 hPa u results. 

ii) Inter-X r48 U 200 

The t4 8 U 200 inter-x cross spectra are shown in Fig. 30b. 
Compared to the t00 diagram (Fig. 30a) , the gross patterns are 
roughly similar, with some notable differences. East of 150°E 
the coherence is higher between 12.5°N-20°N, and lower in the 
equatorial belt and to the south. The eastward propagating 
signals also expand southward from that indicated in the TOO 
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chart, to 12.5°N. This result is similar to the comparison of 
the t48 vs TOO V 200 inter-x spectra. It seems that in the area 
of forecasted eastward propagation the forecasted field has 
better inter-x organization downstream. To the south the 
reverse appears true. Here the waves are propagating westward, 
and upsteam from the base longitude the wave structure is less 
organized than indicated by the TOO analysis. 

Also in Fig. 30b the inter-x cross spectra in the entire 
equatorial zone show a tendency towards westward phase 
propagation but very small phase differences over all 
longitudes. Thus a zonal scale much larger than the wavenumber 
three shown in the TOO U 200 is indicated. This is a poorer 
performance compared to that of u at 850 hPa (Figs. 28a, b), 
where the 48-h forecast successfully reproduces the 120° 
longitude wavelength found in the TOO inter-x cross spectra. 



5. Meridional structure of the synoptic wave disturbances 

In this section the meridional structure for each field is 
examined by comparing the cross spectra between time series at a 
grid point and that at 10°N of the same longitude. Here 10°N is 
chosen to be the base latitude. The results will be referred to 
as the inter-y cross spectra, and are plotted as vectors 
representing the phase differences and coherence squares in the 
same way as the inter-x cross spectra. 



a. Meridional cross spectra of v g^n 

Fig. 31a displays the inter-y cross spectra for TOO vg 5 Q. South 
of 10°N, most vectors rotate clockwise with increasing phase 
difference southward. This phase characteristic is especially 
prominent in the western Pacific east of 140°E. The positive 
phase angle indicates that the v fluctuations at 10°N lags those 
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at latitudes to the south. Because the waves propagate westward 
over this area, a northeast-southwest tilt is indicated. The 
phase difference approaches a quarter-cycle between 10°S and 
10°N. To the north of 10°N most of the vectors in the western 
tropical Pacific show an opposite phase difference, which again 
implies a northeast-southwest tilt. While there are no direct 
meridional cross-spectra in the RAOB data available to verify 
this result, this structure is consistent with all previously 
observed meridional structure of the easterly waves in the 
northern tropics (Chang et al., 1970, 1977; Reed and Recker, 
1971; Burpee, 1972, 1974; Reed et al, 1977). The T48 Vg 5 Q 

inter-y cross spectra are shown in Fig. 31b. They resemble the 
patterns of Fig. 31a very closely. The only difference is that 
the northeast-southwest tilt in the equator-15 °N band is 
somewhat less prominent compared to rOO. 

b. Meridional cross spectra of u g 50 

The rOO ug 5 Q meridional structure (Fig. 32a) also indicate a 
tendency of the northeast-southwest tilt in the tropical western 
Pacific, but the phase difference is less prominent or even 
changes direction south of 10°N. The tilt cannot be traced 
continuously south of the equator, and in the region east of 
150°E the phase is about a half-cycle out of phase with the 
northern tropics. The out-of-phase relationship indicates that 
ug 5 Q is anti-symmetric about the equator, a structure that is 
consistent with the mixed Rossby-gravity wave mode. Coupled 
with the result of a longer zonal wavelength in the equatorial 
belt particularly in the u field, it seems that in the 
equatorial belt the 5-day u fluctuations are mainly 
manifestations of the mixed Rossby-gravity waves. 

The inter-y cross spectra for r48 ug 5 Q are shown in Fig. 32b. 
Compared to TOO, the coherence squares south of 10 °N decrease 
over most of the longitudes except between 150°E-160°E. 
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c. Meridional cross spectra of V 200 



In the TOO V 200 inter-y structure (Fig. 33a) the 
northeast-southwest tilt between the equator-15°N is seen, but 
with a notable exception. Between 130°E-145°E the meridional 
tilt is northwest-southeast. The r48 V 200 results (Fig. 33b) , 
however, show a northeast-southwest tilt throughout the 
equatorial western Pacific. 



d. Meridional cross spectra of U 200 

Fig. 34a shows the northeast-southwest tilt throughout the 
western Pacific in the TOO U 200 inter-y cross spectra. In the 
t48 U 200 results, shown in Fig. 34b, this tilt is much less 
prominent . 

6. Vertical structure of the synoptic wave disturbances 

The vertical structure of the synoptic wave disturbances is 
determined from cross-spectral analyses between different 
vertical levels. Figs. 35-38 show the results at the 5-day band 
(averaged over 4.6-, 5.3- and 6.4-day) with v at 850 hPa as the 

base time series and v at 700, 500, 300 and 200 hPa, 

respectively, as the other series. Each of these figures have 
two panels, panel a is for t00 and panel b is for t48. For 
these cross spectra a positive phase difference indicates that 
v 850 lags the other series. 

Fig. 35a shows that the coherence squares are generally high 
between t00 Vgso and V 700 for the 5-day band, apparently due to 

the closeness of the two levels. North of 20 °N where the waves 

propagate eastward, the phase difference is negative and 
indicates that vg 5 Q leads v-jqq. This means that there is a 
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westward tilt with height, as may be expected for the 
midlatitude baroclinic systems. In the tropical westward 
propagating region most of the phase difference is small, 
indicating a rather small tilt. The main exception is over the 
equatorial region east of 150°E, where the phase differences 
indicate a eastward tilt with height. Fig. 35b shows that t48 
reproduces the TOO structure very well. 

Fig. 36a shows the vertical tilt between Vgso and V 500 at TOO. 
The westward tilt north of 15 °N is more prominent, and the 
eastward tilt is visible in the southern equatorial region east 
of 150 0 E . In addition, west of 140°E and east of the 
Philippines an eastward tilt of up to 0.1 cycle can be found as 
far south as 10 °N and slightly beyond. This is the region of 
Yap and Koror where previous spectral analysis found a westward 
vertical tilt in the easterly waves during summer (Chang et al, 
1970, 1977) . Holton (1970) has proposed in a numerical 
simulation that the vertical shear of the tropical easterly 
waves are related to the vertical shear of the time-mean flow. 
The difference in the vertical tilt between the present study 
and the previous ones may be attributable to the difference in 
the vertical shear, which changes from a easterly shear during 
the northern-summer southwest monsoon to a westerly shear during 
the northern-winter northeast monsoon. The coherence squares 
are below the significant level southwest of 10°N, 140°E and in 
the South China Sea. The lack of coherence between 850 hPa and 
500 hPa in the latter region is expected, as this is the period 
where synoptic timescale fluctuations are heavily influenced by 
winter monsoon surge events whose vertical structure is limited 
to below 700 hPa (Chang et al., 1983; Boyle and Chen, 1987). 

Fig. 36b shows that again the t 48 results simulate the TOO 
structure very well. It may be noticed that between the equator 
and 10°N, the coherence values for T48 are slightly higher than 
TOO to the west of 160°E and slightly lower to the east. 
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Namely, west of 160°E the vertical structure of the 5-day waves 
in the 500-850 hPa layer is somewhat more clearly defined in the 
r 48 forecast field than the rOO analysis, and east of 160°E it 
is somewhat less clearly defined in the r48 field than the TOO 
field. As will be discussed in the following paragraphs, this 
westward shifting of the area of vertical coherence from TOO to 
r48 in the equatorial central-western Pacific becomes more 
prominent when the cross spectra between Vg 5 Q and upper 
tropospheric v are examined. 

Fig. 37a shows the TOO vertical structure between V 350 and V3oo* 
Unlike previous diagrams (Figs. 35-36) that show significant 
vertical coherence over a large longitudinal span in the lower 
troposphere over the equatorial western Pacific, the vertical 
coherence between lower and upper troposphere is restricted to 
just east of 160°E. Between the equator and 10°N, the westward 
propagating 5 -day waves show a westward tilt with height of 
around one quarter cycle phase shift between 850 and 300 hPa. 

As the waves propagate westward across 160°E, signals of 
vertical coupling disappear. Fig. 37b shows the vg 5 Q - V 300 
structure as represented by the t48 forecast fields. A strong 
contrast can be seen between Figs. 37a and 37b. East of 160°E, 
the significant vertical coupling signals analyzed by TOO data 
in the equatorial region (Fig. 37a) weaken somewhat in the t48 
forecasts (Fig. 37b) . On the other hand, west of 160°E the t48 
data generate a significant vertical structure that closely 
couple the two levels. The structure is a westward tilt with a 
phase shift that increases progressively toward west. This 
vertical coupling is not present in the TOO data. 

The difference in the lower-upper troposphere vertical coupling 
over the equatorial Pacific between the TOO analysis and t48 
forecast is seen even more clearly in the vgjg - v 200 cross 
spectra (Fig. 38) . In the equator-10°N band the TOO result 
(Fig. 38a) shows a weak vertical coupling east of 160°E, with 
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some grid points having marginally significant coherence values 
indicating a westward tilt. West of this longitude, between 
140°E-160°E, the two levels are basically not coupled. West of 
140 °E the northern equatorial region shows a limited area of 
high coherence and large out-of-phase tilt between Vg 5 Q and 
v 200* For 748 ( F ig* 38b) this variation with respect to the 
160 0 E longitude is reversed. East of 160°E the weak vertical 
coupling signals in the TOO analysis is now even weaker, with no 
significant coherence on either side of the equator. West of 
160°E where the TOO data show no vertical coupling, r48 shows 
large coherence values and a clear westward tilt with height. 

As is shown in the r48 Vg 5 Q - V 300 results, this tilt also 
increases progressively as the waves propagate toward the west, 
and becomes nearly one half cycle out of phase at 150°E. 

The vertical structure within the upper troposphere can be 
determined from Fig. 39, which shows the vertical cross spectra 
between V 300 ~ v 200' with V 200 as th e t>a se series. Here both 
rOO (Fig. 39a) and t48 (Fig. 39b) show a significant westward 
tilt structure between the two levels in the equatorial 
central-west Pacific. While the two diagrams resemble each 
other reasonably well, it may still be seen that the t48 data 
appear to move the maximum coherence westward from near 170°E in 
the TOO diagram to around 152. 5°E. The very modest vertical 
tilt in the latter area is also replaced by the westward tilt 
seen in the upstream region of the TOO result. 

To summarize. Figs. 37-38 indicate that in the equatorial 
central-western Pacific the t48 forecasts appear to generate 
artificial vertically-organized wave structure that couples the 
lower and upper troposphere downstream from where a weakly 
coupled structure is contained in the TOO analyses. The 
artificial vertical structure gives rise to a progressively 
larger westward tilt with height from 160 °E westward to 150°E 
and beyond. Figs. 36 and 39 indicate that this tendency of 
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downstream propagation of the wave structure by the t 48 forecast 
can also be detected within the lower troposphere (850-500 hPa) 
and within the upper troposphere (300-200 hPa) , respectively. 

The location of the Caroline Island stations along 7°-7.5°N 
allows an evaluation of this result by comparing the rOO and r48 
vertical cross spectra with those from the RAOB ' s to the west 
and to the east of the 160°E longitude. Fig. 40 shows the 
coherence squares between vg 5 Q and v at 700, 500, 400, 300 and 
200 hPa for the frequency of (5.3 day) -1 for RAOB, rOO and t 48, 
respectively, at the stations Majuro, Ponape and Truk. Here 
Majuro (7.1°N, 171. 4°E) is east of the 160°E longitude, and 
Ponape(7°N, 158. 2°E) and Truk (7.5°N, 151. 9°E) are to the west 
of 160°E. At all three stations the vg5o _ V7oo coherence values 
are high, and at each station the r48 and rOO coherence squares 
are nearly identical. For cross spectra between vg 5 Q and v at 
500 hPa and higher levels, the T48-T00 coherence comparisons 
show the same anomalous patterns east and west of 160°E as found 
in Figs. 36-38. At Majuro, the r48 values are either comparable 
to or less than the rOO values, indicating that the t 48 data 
forecast a slightly less degree of organization for the 5-day 
oscillations. The reverse is true at Ponape and Truk, where the 
coherence values indicate, in general, a significantly higher 
degree of deep tropospheric coupling in the r48 data that are 
not present in the rOO analysis. 

At all three stations and levels the rOO coherence square is 
always larger than the corresponding RAOB value, which may be 
expected due to the model's three-dimensional objective 
analysis. The "ground truth", or observed, degree of coupling 
between lower and upper troposphere may be represented by the 
RAOB coherence squares of vg5Q-V2oo and vg5Q - V3oo* These values 
are higher for Majuro than for Ponape and Truk, to the extent 
that for Vg5o - V3oo the coherence is above significance level at 
Majuro but below it at Ponape and Truk. Thus the RAOB confirms 
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the rOO result that the vertical coupling is significant only 
east of 160°E. 

If the vertical structure difference across 160°E between r48 
and TOO is indeed a result of an artificial downstream 
propagation by the forecast model, we may expect that an 
intermediate structure will show in the r24 data. Fig. 41 shows 
the r24 vertical cross spectra between V 350 as the base series 
and v at (a) 300 hPa and (b) 200 hPa, respectively. A 
comparison of Fig. 41a with Figs. 37a, b readily shows that, in 
the region east of 160°E between equator-10°N, the Vg50 -V 300 
coherence values at r24 are higher than those at TOO and lower 
than those at t48. The vertical phase shift for t24, an 
increasingly larger westward tilt with height, is consistent 
with that at t48. The same conclusions can be drawn in the 
comparison of the vertical cross spectra at t24 (Fig. 41b) with 
those at TOO (Fig. 38a) and T48 (Fig. 38b) . 

7. Thermal structure of the synoptic wave disturbances 

The thermal structure of the 5-day disturbances may be examined 
by a cross-spectral analysis between the meridional velocity and 
temperature (T) at various levels. However, because of the 
smallness of the temperature fluctuations for tropical synoptic 
scale motions, the temperature signals are often too weak to 
produce cross spectra that are significant (Chang et al., 1970). 
Due to the uncertainties involved, we will only look at the v-T 
cross spectra at the six grid points that are closest to a RAOB 
station in the Caroline Islands network. Table 2 shows the 
coherence squares and phase differences between V 350 and 
temperature at 925, 850, 700, 500, 400, and 300 hPa. Here only 
those results with a coherence square > 0.14 (95% significance 
level) are formally included, but those with a coherence square 
> 0.11 are also listed. It is seen that only the three eastern 
stations, Majuro, Truk and Ponape, produce substantial numbers 
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of significant results, and they are mainly for the 
upper-tropospheric temperatures. 



Table 2. Coherence square and phase difference between vg 5 Q as 
the base series and temperature at different levels for the 
Caroline Islands region. Values associated with CH 2 > .14 (95% 
significant level) are shown as regular types. Values 
associated with .11 <CH 2 < .14 are shown in italics. Values 
associated with CH 2 < .11 are not shown. 
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Table 2. (continued) 



Koror 

( 134 °E, 7 . 3 °N) 
RAOB TOO r48 
T 3 oo CH 2 - 
Phase - 



Yap 

( 138 °E, 9 . 5 °N) 
RAOB TOO T48 



Guam 

(145°E, 13 . 6 °N) 
RAOB TOO r48 



T 400 CH 2 - 
Phase - 



T 5 00 CH 2 

Phase 



.11 .11 

55 ° 76 0 



T700 CH2 

Phase 



. 17 
60 



. 14 
86 



T 8 50 CH 2 
Phase 



. 62 
63 



.12 

62 ° 



T 925 CH 2 
Phase 



.11 - .20 

- 104 ° - 53° 



At Ponape and Truk both the TOO and t 48 results show an 
approximate one-quarter cycle phase difference between vgso and 
temperature at the 300 and 400 hPa levels. This phase 
difference is somewhat smaller (about 1/6 cycle) with respect to 
T 500 • The positive phase indicates that maximum temperature 
leads maximum V 350 , thus placing the upper-tropospheric warm 
core in the 850 hPa trough region. Therefore the model analysis 
and forecast data both show a "warm core" structure for the 
5-day disturbances, and this agrees with previous results based 
on station data (Chang et al, 1970, 1977). However, the RAOB 
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v 850 -t cross spectra give significant results only at Ponape. 
Here the phase differences indicate a warm core structure around 
500 hPa, but the result for T 200 i s below the significance 
level. The warm core structure would, hydrostatically, suggest 
a weakening of the wave trough with height. Compared to the 
vertical cross spectra between v at different levels, which 
shows very little tilt below 500 hPa and a westward tilt above, 
we may sketch a schematic diagram of the vertical cross section 
of the 5-day disturbances as shown in Fig. 42a. At Majuro t 48 
data do not give significant Vg 5 Q-T results in the middle-upper 
troposphere. The TOO results in the middle levels are 
consistent with Truk and Ponape, but the RAOB data indicate that 
while the disturbances are still warm core, the phase difference 
between Vg 5 Q and middle-upper level temperature is slightly 
smaller. The schematic structure shown in Fig. 42a is basically 
applicable at Majuro, except that the upper warm core is 
slightly shifted to the east and that the upper tropospheric 
westward tilt is slightly stronger according to the vertical 
cross spectra between v at different levels shown earlier (Fig. 
37a) . 

At the three eastern stations, evidences of a model that is 
trying to organize the middle-upper tropospheric wave structure 
downstream in the direction of propagation can again be found. 
Basically, west of 160°E at Majuro the r48 coherence is weaker 
than TOO, while east of 160°E at Truk and Ponape the reverse is 
true. 

At the lower levels the few significant cross spectral results 
give a confusing picture of the thermal structure near the 
surface. For TOO, both Ponape and Truk give a positive 
quarter-cycle phase difference between Vg 5 Q-T 925 , indicating a 
temperature maximum that is approximately in-phase with the 
low-level trough. This warm core at 925 hPa would place the 
maximum amplitude of the wave trough near the surface in order 
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to satisfy hydrostatic balance. On the other hand, the t 48 
results at Majuro, Truk and Yap indicate that the phase of 
v 850~ t 925 about -1/4 cycle, or a cold core at 925 hPa. 
Unfortunately, the corresponding RAOB information is not 
available because 925 hPa is not a standard data level. 

Although previous observational studies by Chang et al (1970) 
suggests that the tropical easterly waves over the western 
Pacific tend to have a cold core structure, mainly during 
summer, near the surface, we will have to view the present 
result as showing an artificial model effect that converts the 
analyzed low-level warm core to cold core in the r48 forecast. 
There are several possible mechanisms that may cause a low-level 
cold core structure. If forced lifting such as CISK is present 
below a cloud base that is too high, adiabatic cooling will 
occur. If precipitation is intense through an unsaturated 
environment near the surface, the re-evaporation of the falling 
rain will cause cooling. Another possibility is that frictional 
damping is too strong and dissipates too much kinetic energy 
near the surface, resulting in a weaker wave amplitude there and 
the hydrostatic constraint will produce a cold-core trough. 

At the other three stations, only a few cross-spectra of vg 5 Q-T 
above the significance level are obtained in the lower 
troposphere, and none in the upper troposphere. At Yap and 
Koror the r00 data show a positive phase of about 1/5 cycle for 
t 700' suggesting that the middle tropospheric warm core over the 
850 hPa trough extends to lower levels at these western 
stations. This seems to be supported by the RAOB result for 
Tg 5 Q at Yap. The V 850 -T 925 result for TOO at Guam, and the less 
significant result for t48 at Yap, are consistent with the 
difference between model analysis and forecast at the three 
eastern stations. Even though the upper tropospheric 
temperature signals are not discernible, a schematic diagram for 
the vertical structure at these western stations is given in 
Fig. 42b for comparison with the other stations. 
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8. Velocity potential structure 



It would be interesting to determine the propagating structure 
of the divergence field associated with the 5-day oscillations. 
However, the divergence itself normally contains substantial 
small-scale, often very noisy, features, so the velocity 
potential is used in the analysis. The velocity potential is 
computed following the procedure discussed by Boyle and Chang 
(1984). Because it is equivalent to a field that is 
approximately proportional to the spatially-smoothed divergence, 
it may only represent scales that are larger than the synoptic 
disturbances. Since the divergent component of the motion is 
usually not well observed, the divergence contained in the 
objective analysis of NWP systems is normally heavily influenced 
by the model forecast through the four-dimensional data 
assimilation (Trenberth, 1992) . Comparison with a non-NWP 
assimilated data set indicates that the assimilated data tends 
to smear out smaller scale features in the velocity potential 
field (Chang et al 1991) . Thus the divergence field associated 
with the tropical synoptic scale disturbances is 

under-represented in the velocity potential. Nevertheless, the 
horizontal structure of the velocity potential is determined 
using the inter-x and inter-y cross spectral analyses mainly to 
compare the 748 forecasts with the TOO analyses within the 5-day 
band, which should be relevant to the evaluation of the forecast 
for the planetary scale that seems to dominate the equatorial 
region. The 150°E and 10°N are again used as the base longitude 
and base latitude, respectively. 

a. Longitudinal cross spectra of velocity potential 

The 850 hPa velocity potential inter-x cross spectra shown in 
Fig. 43a shows a planetary scale eastward propagation pattern 



39 



extending from the northern boundary southward to cover most of 
the northern tropics. Westward propagation, with a very large 
zonal scale, can only be seen at the equator and to the south. 
Thus any divergence field associated with the 5-day synoptic 
easterly wave disturbances observed in the vgso/ which shows 
westward propagation at and south of 15°N (Fig. 27a) , is 
apparently masked by the midlatitude-based eastward propagating 
divergence system north of the equator. The southern equatorial 
westward propagating velocity potential pattern also possesses a 
much larger zonal scale than the wavenumber-three waves observed 
in the ug 5 Q time series there. 

Fig. 43b shows the inter-x structure of the r48 velocity 
potential at 850 hPa. It is a good reproduction of rOO, 
probably more a result of the fact that rOO is substantially 
influenced by the model forecast than a good 48-h forecast of 
the actual field. The difference between Figs. 43a and 43b is 
mainly in the equatorial latitudes, where the coherence over the 
South China Sea is increased. A larger (quarter-cycle) phase 
difference between the western Indonesia region and 150 °E 
indicates that the westward propagating waves have a wavenumber 
two zonal scale, which is shorter than that indicated by the rOO 
analysis in Fig. 43a. 

The inter-x structure of the rOO velocity potential at 200 hPa 
is shown in Fig. 44a. Here the indicated zonal scales are 
comparable to those shown at 850 hPa, but the westward 
propagating pattern is identifiable in the tropical band south 
of 15 °N . Thus at the 5-day frequency band, the region between 
the equator and 15°N is dominated by an eastward propagating 
midlatitude divergence system in the lower troposphere, and by a 
westward propagating tropical divergence system in the upper 
troposphere. Fig. 44b, which is the inter-x structure of T48 
velocity potential at 200 hPa, indicates that the upper 
tropospheric pattern is well forecasted north of the equator 
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only. In the southern equatorial region, the westward 
propagation is poorly represented, particularly toward the east 
of 150°E where the phase vector indicates a slight eastward 
propagation. 

b. Meridional cross spectra cf velocity potential 

The rOO velocity potential inter-y cross spectra at 850 hPa 
(Fig. 45a) demonstrate high coherence values throughout the 
domain, indicating that the 5-day fluctuations are of a very 
large north-south scale. Again, this is characteristic of 
four-dimensional assimilated data where the divergence field is 
significantly affected by the numerical model forecast. Over 
the tropical western Pacific between the equator and 20°N, the 
fluctuations are in phase. Outside of this region, the 
meridional tilt is predominately northeast-southwest. This is 
true in the southern equatorial region also, where the negative 
phase vectors imply a northeast-southwest tilt because the zonal 
propagation is westward (Fig. 43a) . Another view of this is 
simply that there is a phase propagation southward throughout 
the domain. Since this is occurring at the synoptic timescale, 
it may be a manifestation of the southward cold surges from the 
East Asian continent observed in the northeast monsoon region 
(Chang et al, 1983, Chang and Chen, 1992, Chang et al. 1993). 

The fact that the southward propagation can be found south of 
the equator indicates that the cold surge associated convergence 
can penetrate into the southern hemisphere. This general 
structure is reproduced by the r48 850 hPa velocity potential 
(Fig. 45b) . 

The TOO velocity potential at 200 hPa (Fig. 46a) shows a 
meridional structure that is different from that at 850 hPa . In 
the northern eastward-propagation region east of 120°E, the tilt 
is predominantly northwest-southeast, which is opposite to the 
850 hPa structure. South of 10°N in the westward-propagation 
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region east of 120°E, the tilt is northeast-southwest, again 
opposite to the 850 hPa structure. The two regions combine to 
show a northward propagation throughout the domain. Previous 
studies of winter monsoon cold surges (Chang and Lau, 1982, 

Chang and Chen, 1992) reported enhancement of upper tropospheric 
northward divergent wind following each cold surge event. The 
present result indicates that the divergence field itself is 
propagating northward at the synoptic timescale. 

The meridional tilts or propagation almost disappear entirely in 
the r48 time series (Fig. 46b) . This difference between the t48 
and rOO 200 hPa velocity potential is somewhat surprising, since 
one would expect that the four-dimensional assimilation process 
would strongly influence the divergence field of the model 
forecast. It appears that in this case the analyzed rOO 
divergence field contains observed information in the meridional 
structure that is independent from the model forecast, and that 
the forecast divergence is unable to reproduce this information. 



9. A proposed scheme to evaluate systematic errors in model 
analysis and forecast of tropical disturbances 

The previous sections show that time series analysis may be used 
to evaluate a model's ability to represent tropical disturbances 
during a season. In order to simplify the evaluations, we 
propose a scheme to use spectral characteristics to measure the 
systematic errors of a NWP model with respect to the climatology 
of the tropical oscillations within a given season. The scheme 
includes three error parameters, all for a pre-selected 
frequency band. The model output climatology may be the 
analysis or the forecast, and the verification climatology may 
be the radiosonde observations or the analysis. 

1) Coherence square (Coh2): This is simply the coherence square 

between the model output time series (analysis or forecast) and 
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the verification time series. The higher the Coh2 , the more 
correlation between model output and verification. 

2) Variance difference (VDif) : This measures the errors in 

estimating the amplitude of the oscillations: 

VDif = 2n x of x 2^ [(SVm) 3 ? - (SVy)^] / n, 

Here £>f is the frequency interval in the spectral analysis, £V m 
is the summation of variance (power) at each frequency point 
within the frequency band for the model output time series, £V v 
is the summation of variance at each frequency point for the 
verification time series, and n is the number of frequency 
points within the frequency band. A positive VDIF indicates 
that the model output overestimates, and vice versa. A smaller 
VDIF indicates a more accurate amplitude estimate of the 
oscillations . 

3) Phase shift (PS): This measures the phase shift between the 

model output and the verification: 

PS = x P 

Here cup is the phase difference in cycles (°/360°), and P is the 
averaged period within the frequency band. Thus PS gives the 
phase difference in time units. A smaller PS is more accurate. 

The scores for u and v components for the 5-day band at 200 hPa 
and 850 hPa are computed at each grid point closest to a 
verifying radiosonde station. Figs. 47-50 show the results from 
TOO to r 72 at 24 h intervals for each of the stations in the 
Caroline Island region, with the averages computed using root 
mean squares. In general, the coherence squares show clearly a 
degradation with time. At 200 hPa the average CH2 decreases 
from about 0. 6-0.7 at r00 to less than 0.2 for both u (Fig. 47) 
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and v (Fig. 48), with the largest decline occuring during the 
first 24 hrs. At 850 hPa the values are generally higher, with 
the average CH2 declining from 0.8 to near 0.5 for both u and v 
(Figs. 49-50, respectively) , indicating that the model analysis 
and forecast retain much more the observed synoptic oscillation 
climatolgy in the lower troposphere than the upper troposphere. 

Because root mean squares are used, the average variance 
difference and phase shift are both positive. (Yap is excluded 
from the average VDif due to the large discrepancies between rOO 
and RAOB vgsoO Individual station VDif scores show that at 200 
hPa there is a tendency of decreasing variance in the model 
outputs such that by r 72 the u and v model climatology has an 
amplitude about 1.5 m/s less than observed. The change at 850 
hPa is less systematic. The phase shift varies greatly among 
the stations. At 200 hPa the u results (Fig. 47) show that the 
two western stations (Koror and Yap) have a positive phase error 
at t 48 and r 72 of around 0.5 day, while the four eastern 
stations show negative phase errors average about -0.5 day. 

Other diagrams have less systematic patterns. 

10. Summary and concluding remarks 

In this study we use spectral analysis to diagnose the structure 
of synoptic timescale wave disturbances over the tropical 
western North Pacific during the 1991-1992 winter with three 
data sets. The data are the radiosonde observations in the 
vicinity of the Caroline Islands, and the initial (r00) analyzed 
and the 48-h (t4 8) forecast fields from the Navy Operational 
Global Atmospheric Prediction System (NOGAPS) . The primary 
purpose is to evaluate the ability of the model analysis and 
forecast to represent the characteristics of the tropical 
synoptic disturbances within a frequency band centered at (5 
day) -1 . This is done through a comparison of the wave 
climatologies produced by the rOO and r48 model data with that 
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determined from radiosonde data. The comparison of the r 00 and 
r 48 model climatologies also allows a study of the model's 
systematic errors in representing tropical synoptic 
disturbances . 

As is the case with many previous time series analyses of 
tropical synoptic disturbances, oscillations in the meridional 
wind component, especially at 850 hPa , are used as the main 
indicator of disturbance activity. In many cases two spectral 
power "peaks" are found in the vg 5 Q radiosonde data: a 
low-frequency concentration of power encompassing the seasonal 
and intra-seasonal oscillations, and a synoptic peak centered 
around five days that is the focus of this study. While the TOO 
data, in general, resemble reasonably well the spectral 
distribution of the radiosonde results, the t 48 forecast often 
has markedly larger discrepancies. On the other hand, 
cross-spectral analysis of the meridional wind time series in 
the horizontal domain demonstrates good agreement on the general 
structure of the synoptic waves among the TOO analysis, the t 48 
forecast and the radiosonde observations. As examples: 

1) Based on the time series of meridional wind component at 850 
hPa, the synoptic easterly (westward propagating) waves are well 
defined in the 135°E-180° region in both the radiosonde and the 
model data. The wavelength is » 33° along 10°N and somewhat 
longer (38°) along the equator. These observed characteristics 
are very well represented by both the TOO analysis and the t 48 
forecast . 

2) For regions south of 20°N in the western Pacific, the 
easterly waves exhibit a northeast-southwest tilt in the t 00 
v 8 5 o field, a feature frequently observed by previous studies 
based on radiosonde data. This is well reproduced by the t48 
Vg 5 Q time series. 
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3) The rOO v at 200 hPa shows an easterly propagating wavelength 
of 50-60° longitude., which is reasonably reproduced by the r48 
V 200 data . 

The horizontal cross spectra in both the u component and the 
velocity potential are different from the v component results. 
The synoptic wave structure is not clear in the u time series, 
which has longer zonal scales. The rOO data show a 
westward-propagating wavelength near wavenumber three for the 
5-day oscillations in the equatorial zone at both 850 and 200 
hPa. This wavelength is reproduced by r48 u only at 850 hPa . 

At 200 hPa, the r48 u gives an almost in-phase relationship 
across all longitudes in the equatorial belt. Similarly, at the 
5-day time scale, the analyzed velocity potential field is of a 
larger spatial scale than the synoptic scale shown by the v 
data. The implied large-scale divergence propagates southward 
at 850 hPa and northward at 200 hPa , both of which appear to be 
manifestations of the winter monsoon cold surges observed by 
Chang and Chen (1992). The t48 forecast can reproduce only the 
southward propagation at 850 hPa. It is unable to reproduce the 
northward propagation of the upper-level divergent flow. Since 
the synoptic timescale upper-level divergent return flow has 
been found to be related to the acceleration of the East Asian 
jet streak (Chang and Lau, 1982; Chang and Lum, 1984), this 
inability may cause the model to underforecast the strengthening 
of midlatitude baroclinic systems. 

The vertical structure of the 5-day wave disturbances in the 
Caroline Islands determined from the r00 data is in good 
agreement with that determined from the radiosonde reports. In 
the eastern part of the station network the disturbances have 
very little vertical tilt in the lower troposphere, in agreement 
with the structure observed there previously by Wallace and 
Chang (1969). In the upper troposphere they show a westward 
tilt with height, which is in agreement with the 
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upper-tropospheric structure observed in the central Pacific 
during the northern summer by Chang et al. (1970) , although 
these central Pacific waves have a strong eastward tilt with 
height in the lower troposphere. In the western part of the 
station network near Yap and Koror there is a eastward tilt with 
height from lower to upper troposphere. This is opposite to the 
vertical tilt observed by Chang et al. (1970) for these stations 
during summer of 1964. However, the effect of the vertical 
shear of the time-mean zonal wind (Holton, 1970; Chang and 
Miller, 1977; and Shapiro et al, 1988) can be used to interpret 
this difference consistently as a result of the difference 
between the winter-northeast and summer-southwest monsoons. 

This difference is most pronounced in the lower troposphere, and 
over the western part of the western Pacific. 

West of 160°E between the equator and 10°N, the r48 times series 
tend to produce too strong a wave structure downstream from a 
region where vertically-organized oscillation is observed. This 
may be due to the equatorial dynamics contained within the model 
that in some way propagates the observed wave structure 
downstream into a region where it is does not exist. It is 
fortunate that we have reliable radiosonde stations in the 
Caroline Island region to verify the analysis. Namely, we can 
ascertain that 1) the forecast climatology is artificial, and 2) 
the analysis climatology incorporates observed data to correct 
the model forecast first guess. In fact, the TOO to t 48 
westward shift of the maximum coherence between V 200 ” v 300 i s 
approximately 17.5° longitude (Fig. 39). This is only slightly 
slower than the phase propagation speed indicated by the V 200 
zonal cross spectra, which is 50° longitude per five days, or 
about 20° per 48 h. For the lower levels the westward shift is 
about one-third less, corresponding to the fact that the 
propagation speed at 850 hPa is slower, at 33° per five days. 
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It is possible that this problem of artificial generation or 
propagation of equatorial waves by the numerical models may 
occur elsewhere in the tropics and during other seasons. 

Because of lack of conventional data, many tropical studies use 
four-dimensionally assimilated data that are analyzed by 
operational NWP models. This use of model analysis as real data 
is especially risky over a data void region, since any 
model-generated equatorial wave structure will most likely be 
present in the analysis data as they cannot be corrected by real 
observations . 

The 5-day wave is shown to be warm core above the 850 hPa trough 
in the middle upper troposphere. Therefore it may be 
self-sustaining energetically if synoptic-scale rising motion is 
in phase with the low-level trough. However, the rOO analysis 
and the r48 forecast climatologies give different thermal 
structures near the surface, with the former showing a warm core 
and the latter showing a cold core. This rather dramatic 
difference is perhaps due to the model's cumulus and boundary 
layer parameterizations that produce a cold core near the 
surface. Several model mechanisms that may be involved are 
suggested. They include a higher level of cloud base formed by 
a large-scale process and thus some artificial adiabatic cooling 
below, over-intense rainfalls that re-evaporate in the boundary 
layer, or too strong a lower boundary dissipation effect. 

The zonal velocity component and velocity potential both show a 
much larger wavelength for the 5-day oscillations than that 
deduced from Vg 5 Q, so the synoptic spatial-scale disturbances 
are masked by the planetary scales in these fields. For TOO, 
the velocity potential over the northern tropical western 
Pacific and the Southeast Aisa shows a southward propagation at 
850 hPa and a northward propagation at 200 hPa . The low-level 
propagation can be traced to the southern tropics. This pattern 
appears to be related to the acceleration of the East Asian 
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local Hadley cell during winter monsoon cold surges suggested by 
Chang and Lau (1982). In the r48 forecast climatology, the 
low-level propagation is represented, but the upper level 
poleward propagation is missing. 

Based on the study of the spectral and cross-spectral 
characteristics, a scheme to quantify the systematic errors in 
tropical oscillations in the NWP analysis and forecast is 
proposed. The scheme uses three scores to measure the errors: 
the coherence square, the variance difference, and the phase 
shift, between the oscillation climatology of the model output 
(analysis or forecast) and that of the verification (radiosonde 
observation or analysis) . For the period of study the loss of 
coherence with time is most visible, especially at 200 hPa where 
the coherence square with radiosonde observations decreases from 
around 0.7 at TOO to less than 0.2 at r72. At 850 hPa the model 
climatology has considerably higher coherence square with the 
observations, decreasing from about 0.8 at TOO to 0.5 at r72. 
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Figure Captions 



Fig. 1. Comparison of vgso spectra at Majuro: a) relative power 
spectra for RAOB (thick) , NOGAPS rOO (thin) and NOGAPS T 48 
(dashed) , b) phase difference (line) and coherence square 
(sqaure) between RAOB and TOO, c) phase difference (line) and 
coherence square (sqaure) between RAOB and r48. 



Fig. 2. 


Same as Fig. 1 except for Ponape. 


Fig. 3. 


Same as Fig. 1 except for Truk. 


Fig. 4. 


Same as Fig. 1 except for Guam 


Fig. 5. 


Same as Fig. 1 except for Yap. 


Fig. 6. 


Same as Fig. 1 except for Koror. 


Fig. 7. 


Comparison of U 850 spectra at Majuro: a) relative power 



spectra for RAOB (thick) , rOO (thin) and t48 (dashed) , b) 
phase difference (line) and coherence square (sqaure) between 
RAOB and rOO, c) phase difference (line) and coherence square 
(sqaure) between RAOB and r48. 



Fig. 8. 


Same as Fig. 7 except for Ponape. 


Fig. 9. 


Same as Fig. 7 except for Truk. 


Fig. 10. 


Same as Fig. 7 except for Guam 


Fig. 11. 


Same as Fig. 7 except for Yap. 


Fig. 12. 


Same as Fig. 7 except for Koror. 


Fig. 13. 


Comparison of V 200 spectra at Majuro: a) relative 



power spectra for RAOB (thick) , r 00 (thin) and r48 (dashed) , 
b) phase difference (line) and coherence square (sqaure) 
between RAOB and TOO, c) phase difference (line) and 
coherence square (sqaure) between RAOB and t48. 



Fig. 14. 


Same as Fig. 13 except for Ponape. 


Fig. 15. 


Same as Fig. 13 except for Truk. 


Fig. 16. 


Same as Fig. 13 except for Guam 
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Fig. 17. Same as Fig. 13 except for Yap. 
Fig. 18. Same as Fig. 13 except for Koror. 



Fig. 19. Comparison of U 200 spectra at Majuro: a) relative power 
spectra for RAOB (thick) , TOO (thin) and t48 (dashed) , b) 
phase difference (line) and coherence square (sqaure) between 
RAOB and TOO, c) phase difference (line) and coherence square 
(sqaure) between RAOB and t48. 



Same as Fig. 19 except for Ponape. 

Same as Fig. 19 except for Truk. 

Same as Fig. 19 except for Guam 

Same as Fig. 19 except for Yap. 

Same as Fig. 19 except for Koror. 

Phase difference (line) and coherence square (square) 



Ponape (P) , Guam (G) , Yap (Y) and Koror (K) for: a) 4.6-day, 
b) 5.3-day, and c) 6.4-day periods. 



Fig. 


20 . 


Fig. 


21 . 


Fig. 


22 . 


Fig. 


23. 


Fig. 


24 . 


Fig. 


25. 


of 


TOO 



Fig. 26. Same as Fig. 25 except for r48 V 350 . 



Fig. 27. Inter-longitude cross spectra for the synoptic band 
(averaged period 5.3 days) for Vg 5 Q with the 150°E data as 
the base series longitude. The phase difference and 
coherence square between each grid point and the base series 
at the same latitude are plotted as a vector, with the 
direction indicating phase difference (northward pointing = 
0 °, increasing clockwise), and the magnitude indicating 
coherence square: a) TOO, and b) r48. The location of the 

Caroline Island radiosonde stations are indicated with the 
first letter of the station name. 

Fig. 28. Inter-longitude cross spectra, same as Fig. 27 except 
for u 850 . 



Fig. 29. Inter-longitude cross spectra, same as Fig. 27 except 
for v 20 o • 
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Fig. 30. Inter-longitude cross spectra, same as Fig. 27 except 
for u 2 oo • 

Fig. 31. Inter-latitude cross spectra for the synoptic band 

(averaged period 5.3 days) for vg 5 o with the 10°N data as the 
base series latitude. The phase difference and coherence 
square between each grid point and the base series at the 
same longitude are plotted as a vector, with the direction 
indicating phase difference (northward pointing = 0°, 
increasing clockwise) , and the magnitude indicating coherence 
square: a) TOO, and b) r48. 



Fig. 32. Inter-latitude cross spectra, same as Fig. 31 except 
for u 850 . 

Fig. 33. Inter-latitude cross spectra, same as Fig. 31 except 
for v 2 oo • 



Fig. 34. Inter-latitude cross spectra, same as Fig. 31 except 
for u 200 . 



Fig. 35. Inter-level cross spectra for the synoptic band 

(averaged period 5.3 days) between v 8 5 q (base series) - V 700 
The phase difference and coherence square between the two 
series at each grid point are plotted as a vector, with the 
direction indicating phase difference (northward pointing = 
0°, increasing clockwise), and the magnitude indicating 
coherence square: a) rOO, and b) r48. 



Fig. 36. Inter-level cross spectra, same as Fig. 35 except for 
v 850 “ v 500 • 

Fig. 37. Inter-level cross spectra, same as Fig. 35 except for 
v 850 ~ v 300 • 

Fig. 38. Inter-level cross spectra, same as Fig. 35 except for 
v 850 “ v 2 00 • 

Fig. 39. Inter-level cross spectra, same as Fig. 35 except for 
V200 (base series) - v 30 o- 
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Fig. 40. Comparisons between RAOB, TOO and r48 inter-level 
coherence squares (vg 5 Q as the base series vs. v at other 
levels) for the 5.3-day synoptic band. See text for details. 

Fig. 41. Inter-level cross spectra, same as Fig. 35 except for 
T24: a) V 350 - V 300 / and b ) v 850 “ v 200* 

Fig. 42. Schematic diagram for the vertical structure of the 
5.3-day synoptic wave at a) eastern stations: Majuro, Ponape 
and Truk; and b) western stations: Guam, Yap and Koror. 
R-ridge axis, T-trough axis, W-warm, C-cold. 

Fig. 43. Inter-longitude cross spectra, same as Fig. 27 except 
for the velocity potential at 850 hPa. 

Fig. 44. Inter-longitude cross spectra, same as Fig. 27 except 
for the velocity potential at 200 hPa. 

Fig. 45. Inter-latitude cross spectra, same as Fig. 31 except 
for the velocity potential at 850 hPa. 

Fig. 46. Inter-latitude cross spectra, same as Fig. 31 except 
for the velocity potential at 200 hPa. 

Fig. 47. The coherence square, variance difference and phase 
shift error scores band between RAOB and model outputs (r00, 
r24, r4 8 , t 7 2 ) for the 5.3-day synoptic band U 200 at grid 
points closest to each of the radiosonde stations. The 
averages of these points are also plotted. 

Fig. 48. Same as Fig. 47 except for V200* 

Fig. 49. Same as Fig. 47 except for ug 5 Q/ and the averages do 
not include Yap. 

Fig. 50. Same as Fig. 47 except for V 350 , and the averages do 
not include Yap. 
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Fig. 1. Comparison of vq^q spectra at Majuro: a) relative power spectra for 
RAOB (thick), NOGAPS tOO (thin) and NOGAPS t 48 (dashed), b) phase 
difference (line) and coherence square (sqaure) between RAOB and tOO, c) 
phase difference (line) and coherence square (sqaure) between RAOB and 
t48 . 
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Fig. 2 
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. Same as Fig. 1 except for Ponape. 
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. Same as Fig. 1 except for Truk. 
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. Same as Fig. 1 except for Guam 



60 



coherence coherence 



YAP ( 9.5N.138.1E) 

PARAMETER = 850 M8 MERIDIONAL WIND 
91111200 - 92031812 
COVERAGE = 90% 




b) 




c) 




Fig. 5. Same as Fig. 1 except for Yap. 
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. Same as Fig. 1 except for Koror. 
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Fig. 7. Comparison of u 850 spectra at Majuro: a) relative power spectra for 
RAOB (thick), rOO (thin) and r48 (dashed), b) phase difference (line) and 
coherence square (sqaure) between RAOB and tOO, c) phase difference 
(line) and coherence square (sqaure) between RAOB and t48. 
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Fig. 8 
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. Same as Fig. 7 except for Ponape. 
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Fig. 9 
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. Same as Fig. 7 except for Truk. 
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Fig. 10. Same as Fig. 7 except for Guam 
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Fig. 11. Same as Fig. 7 except for Yap. 
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Fig. 12. Same as Fig. 7 except for Koror. 
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Fig. 13. Comparison of V 200 spectra at Majuro: a) relative power spectra 

for RAOB (thick), rOO (thin) and t 48 (dashed), b) phase difference (line) 
and coherence square (sqaure) between RAOB and rOO, c) phase difference 
(line) and coherence square (sqaure) between RAOB and t48. 
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Fig. 14. Same as Fig. 13 except for Ponape. 
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Fig. 15. Same as Fig. 13 except for Truk. 
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Fig. 16. Same as Fig. 13 except for Guam 
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Fig. 17. Same as Fig. 13 except for Yap. 
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Fig. 18. Same as Fig. 13 except for Koror. 
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Fig. 19. Comparison of U 200 spectra at Majuro: a) relative power spectra for 
RAOB (thick), rOO (thin) and t 48 (dashed), b) phase difference (line) and 
coherence square (sqaure) between RAOB and rOO, c) phase difference 
(line) and coherence square (sqaure) between RAOB and r48. 
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Fig. 20. Same as Fig. 19 except for Ponape. 



76 



coherence coherence 



TRUK ( 7.5N.151.9E) 

PARAMETER = 200 MB ZONAL WIND 
91111200 - 92031 81 2 
COVERAGE = 86% 



1 .0 

g O.B 

CD 

0.6 

t. 

Q- 0.4 

1 0.2 




b) 




c) 




Fig. 21. Same as Fig. 19 except for Truk. 
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Fig. 22. Same as Fig. 19 except for Guam 
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Fig. 23. Same as Fig. 19 except for Yap. 
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Fig. 24. Same as Fig. 19 except for Koror. 
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Fig. 25. Phase difference (line) and coherence square (square) of r00 V 35 Q 
between Truk as the base series and Majuro (M) , Ponape (P), Guam (G), Yap 
(Y) and Koror (K) for: a) 4.6-day, b) 5.3-day, and c) 6.4-day periods. 
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Fig. 26. Same as Fig. 25 except for t 48 vggg. 



82 



COHERENCE COHERENCE COHERENCE 



(a) 

30N 



20N 



ION 



EQ 



-ICS 



850 mb v tau 0 

COHERENCE/PHASE period = 5.3 days 




100E 120E 140E 160E 180E 

- . 850 mb v tau 48 

COHERENCE/PHASE period = 5.3 days 




1C0E 120E 1 40E 160E 1B0E 



Fig 



27. Inter-longitude cross spectra for the synoptic band (averaged 
period 5.3 days) for vggQ with the 150°E data as the base series 
longitude. The phase difference and coherence square between each grid 



point and the base series at the same latitude are plotted as a vector, 
with the direction indicating phase difference (northward pointing = 0° 



increasing clockwise), and the magnitude indicating coherence square: a) 
tOO, and b) t 48. The location of the Caroline Island radiosonde stations 
are indicated with the first letter of the station name. 
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Fig. 28 



Inter-longitude cross spectra, same as Fig. 27 except for uggg- 
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Fig. 29. Inter-longitude cross spectra, same as Fig. 27 except for V200- 
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Fig. 30. Inter-longitude cross spectra, same as Fig. 27 except for U2QO* 
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. 31. Inter-latitude cross spectra for the synoptic band (averaged 
period 5.3 days) for v 850 with the 10°N data as the base series latitude. 
The phase difference and coherence square between each grid point and the 
base series at the same longitude are plotted as a vector, with the 
direction indicating phase difference (northward pointing - 0 , 
increasing clockwise), and the magnitude indicating coherence square: a) 

TOO, and b) t48. 
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Fig. 32. Inter-latitude cross spectra, same as Fig. 31 except for uqjq. 
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Fig. 33. Inter-latitude cross spectra, same as Fig. 31 except for v 200* 
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Fig. 34. Inter-latitude cross spectra, same as Fig. 31 except for U2QO* 
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, 35. Inter-level cross spectra for the synoptic band (averaged period 
5.3 days) between vggg (base series) - v^qo* P^ ase difference and 

coherence square between the two series at each grid point are plotted 
a vector, with the direction indicating phase difference (northward 
pointing — 0°, increasing clockwise), and the magnitude indicating 
coherence square: a) r00, and b) t48. 
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Fig. 36. Inter-level cross spectra, same as Fig. 35 except for vggQ - v 500* 
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Fig. 37. Inter-level cross spectra, same as Fig. 35 except for vg5Q ~ v 300‘ 
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Fig. 38. Inter-level cross spectra, same as Fig. 35 except for vggQ - v2qo* 
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Fig. 39. Inter-level cross spectra, same as Fig. 35 except for V 200 ( t)ase 
series) - v 300 . 
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Fig. 41. Inter-level cross spectra, same as Fig. 35 except for t 24: a) vggg 
" v 300' and b > v 850 “ v 200 • 



97 



P(mb) 




R T R 

P(mb) 




R T R 



. 42. Schematic diagram for the vertical structure of the 5.3-day 
synoptic wave at a) eastern stations: Majuro, Ponape and Truk; and b) 
western stations: Guam, Yap and Koror. R-ridge axis, T-trough axis, 
W-warm, C-cold. 
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Fig. 43. Inter-longitude cross spectra, same as Fig. 27 except for the 
velocity potential at 850 hPa. 
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Fig. 44. Inter-longitude cross spectra, same as Fig. 27 except for the 
velocity potential at 200 hPa. 
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Fig. 45. Inter-latitude cross spectra, same as Fig. 31 except for the 
velocity potential at 850 hPa. 
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Fig. 46. Inter-latitude cross spectra, same as Fig. 31 except for the 
velocity potential at 200 hPa. 
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Fig* 47* The coherence square, variance difference and phase shift error 
scores band between RAOB and model outputs (r00, t24, t48, t 72) for the 
5.3-day synoptic band u 2 qo at 9 rid points closest to each of the 
radiosonde stations. The averages of these points are also plotted. 
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Same as Fig. 47 except for V200* 
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Fig. 49. 
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Same as Fig. 47 except for UQ 5 Q, and the averages do not include 
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Fig. 50. Same as Fig. 47 except for vggQ, and the averages do not include 
Yap. 
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